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I. INTRODUCTION 

The issue of gravitational stability, how different as¬ 
semblies of massive particles and objects can attain a 
gravitationally-stable state, has been one of the key prob¬ 
lems of astrophysical and cosmological research since the 
late 17th century. Shortly after the publication of his 
monumental Principia , Isaac Newton was led to consider 
the possibility of a spatially infinite universe due to the 
instability of a finite self-gravitating sphere of matter [l]. 
Barely a year after Einstein published his landmark paper 
on the general theory of relativity, he examined the grav¬ 
itational stability of a static, spherically-symmetric uni¬ 
verse [21. That the solution was unstable to perturbations 
led Einstein to include the so-called cosmological term, 
which remains a viable explanation to the accelerating 
recession of far-away Type la supernovae 00, although 
other explanations based on slowly evolving fields are also 
consistent with current data [5]. Moving from the uni¬ 
verse to the stability of compact astrophysical objects- 
the focus of the present manuscript-much depends on the 
object’s specific material composition. In general, the ob¬ 
ject is modeled with an effective equation of state which 
attempts to describe the essential physical processes tak¬ 
ing place in its interior. 

Without presuming to offer here a complete history 
of how gravitational stability of stellar objects has been 
examined over the past decades (for reviews see Refs. 
00 ), we note that an essential aspect of such stability 
is that the object’s binding energy E 5 must be negative- 
definite, Eb = M — Qm < 0. (We take c = h = 1 
throughout, unless explicitily shown.) If E^ > 0 the star 
is unstable to fission. Here, Q is the conserved number of 
particles in the object, either baryon number in ordinary 
stars or the net number of spin -0 bosons in boson stars 
[ 8 ], M is the object’s mass, and m is the mass of the 
particle. 
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It is also well-known that the negativity of the bind¬ 
ing energy is a necessary but not sufficient condition for 
stability: configurations that have E^ < 0 may still be 
unstable under radial perturbations. In order to establish 
the stability of stellar configurations one needs to apply 
perturbations to the effective equations describing the 
self-gravitating matter. Leaving details aside (the inter¬ 
ested reader can consult Refs. 00), the key point is that 
unstable stellar configurations will have exponentially- 
growing radial perturbations characterized by imaginary 
eigenvalues of the perturbed linearized equations. Stabil¬ 
ity conditions are established imposing that the pertur¬ 
bations conserve particle number, as in Chandrasekhar’s 
variational formalism [9]. 

In the present work we examine the stability of self- 
gravitating objects from a very different perspective. In¬ 
stead of the usual perturbative approach, we will apply a 
recently-proposed measure of shape complexity known as 
configurational entropy (CE) [lO] to stellar-like objects, 
investigating the stability of both Newtonian and fully 
relativistic objects. In particular, we will expand the 
results of Ref. mi and apply our formalism to three 
classes of objects: general Newtonian polytropes that 
model non-relativistic and ultra-relativistic white dwarfs 
00 ; neutron stars modeled by an Oppenheimer-Volkoff 
equation of state [ 12 ]; and to boson stars, self-gravitating 
configurations made of charge-conserving complex scalar 
fields [13IIT8] . 

We will show that the configurational entropy, an 
extension of Shannon’s information entropy m to 
spatially-localized mathematical functions based on their 
Fourier transforms, can provide reliable bounds on the 
stability of self-gravitating objects with accuracies of a 
few percent or better. As such, the information-entropic 
method used here provides an alternative approach to the 
study of gravitational stability with broad applicability. 

This paper is organized as follows. In Section II we 
briefly review the formalism and obtain the equations de¬ 
scribing general relativistic, sphericaly-symmetric com¬ 
pact objects. In Section III we review the notion of con¬ 
figurational entropy (CE). In Section IV we apply the 
formalism to cold white-dwarfs, showing how the CE can 
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be used to obtain an estimate of Chandrasekhar’s criti¬ 
cal stability mass for these objects m- We show that 
the mass for polytropes scales inversely with their CE, 
allowing us to relate the mass instability region - a sad¬ 
dle ridge in configuration space - to an equivalent in¬ 
stability region in the object’s configurational entropy. 
In Section V we investigate neutron stars modeled with 
an Oppenheimer-Volkoff (OV) equation of state, showing 
how the CE gives a bound on the compact object’s stabil¬ 
ity consistent with the traditional perturbation method. 
We also show that an inverse scaling relation similar to 
that found for Newtonian polytropes relating the object’s 
mass and its CE holds for OV neutron stars. In Sec¬ 
tion VI we apply the formalism to boson stars with self- 
coupled scalar fields. We show how the CE again gives 
reliable bounds on the compact object’s critical stability 
mass and how the same scaling between mass and CE 
found for polytropes and OV neutron stars is also appli¬ 
cable for these objects. We conclude in Section VII with 
final remarks and a discussion of future projects. 


II. GENERAL FORMALISM 


We consider static, spherically-symmetric configura¬ 
tions with spacetime metric (we follow the conventions 
of ref. [6]), 

ds 2 = B(r)dt 2 — A(r)dr 2 — r 2 (d0 2 + sin 2 6dq i 2 ), (1) 

and take c — h= 1. Einstein’s field equations are 

G„u = -8ttGT mi/ , (2) 

where is the energy-momentum tensor. For Newto¬ 
nian polytropes and neutron stars, we will model stellar 
matter as a perfect fluid with energy-momentum tensor 

Tui; = p(r)g V + \p{r) + p{r)\ U^U„, (3) 

where p(r) is the pressure, p{r) is the energy density 
and U M is the velocity four-vector. Taking the star 
to be at rest, U M has only one non-zero component, 
Uq = —y/B(r). For boson stars, the energy-momentum 
tensor is computed from a Lagrangian density to be de¬ 
fined later. We use the energy density to define the mass 
of the object as 


f R 

M = 4tt p(r)r 2 dr , (4) 

Jo 

where the upper limit of integration, R, is to either 
the object’s radius R, where p(R) = 0 or, for bo¬ 
son stars, to R = oo, although most of the star’s 
mass is concentrated within an effective radius R e ^ = 
J 0 °° p(r)r 3 dr/ f^° p(r)r 2 dr. 

With these definitions, Einstein’s equations can be 


written as: 


if—i — 

A Ar r 2 + r 2 

if— — 

A Br + r 2 r 2 

Bf_ _ 2 p' 

B P + p' 


= 87r Gp; 


= Sir Gp\ 


( 5 ) 


where a prime denotes derivative with respect to the 
radial direction. The last expression is the equation 
for hydrostatic equilibrium. These equations, together 
with an equation of state p(r) = p[p(r)\, are used to 
study a large variety of self-gravitating objects, assum¬ 
ing that A(0) = 1 and B(r oo) = 1. The equation 
involving A{r) and p(r) may be integrated as A(r) = 
[1 — 2GA4(r)/r] _1 , where the mass density function is 
given by A4(r) = J Q r Attt' 2 p(r')dr'. 

As is well-known, Eqs. [5] can be rearranged and, us¬ 
ing the above expression for A(r), the gravitational fields 
A(r) and B{r) can be eliminated to obtain [6 


— r 2 p'(r) = GM{r)p(r) 


47r r 3 p(r 


1 + 


P(r) 

p(r)_ 


M{r) 


1 - 


2 GM(r) 


(6) 


This equation describes self-gravitating stellar config¬ 
urations with general-relativistic corrections in the last 
three terms. We are interested here only in isentropic 
stars, that is, those with a constant entropy per particle 
across the star. Such configurations model very low 
temperature white dwarfs and neutron stars, as well 
as boson stars, which are self-gravitating spin-0 boson- 
condensates. Next we review the main ideas behind the 
configurational entropy measure of spatial complexity 
before we use it to establish stability bounds for all three 
types of configurations. 


III. CONFIGURATIONAL ENTROPY 

Since we are interested in self-gravitating configura¬ 
tions with spatially-localized energy, consider the set 
of square-integrable bounded functions f(x) E T 2 (R) 
and their Fourier transforms F( k). Plancherel’s theorem 
states that 


/ OO /‘OO 

|/(x)|Vx= / \F(k)\ 2 d d k. (7) 

-oo J — oo 

Now define the modal fraction /(k) [TO] . 

_ ijw (8) 

/(k) “ /|F(k)| 2 ^fc’ (8) 

where the integration is over all k where F( k) is well- 
defined and d is the number of spatial dimensions. /(k) 
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measures the relative weight of a given mode k. This can 
also be seen by noting that |F(k )| 2 is proportional to the 
Fourier transform of the two-point correlation function of 
the function /(x), while \F(k)\ 2 d d k is the integrated 
power [20]. For periodic functions where a Fourier series 
is defined, f(k) -» f n = \A n \ 2 / J2\A n \ 2 , where A n is the 
coefficient of the n-th Fourier mode. 

We define the configurational entropy Sc[f] as [W\ 

Self] = - E fn In (fn), (9) 

in analogy with Shannon’s information entropy, Ss = 
OS]- Note that if all N modes k carry 
the same weight f n = 1/7V, the discrete configurational 
entropy has a maximum at Sc = In AT. If only one mode 
is present, Sc = 0. These limits motivate the definition 
of Eq. [9] 

For general, non-periodic functions in the continuous 
interval (a, 6 ), the case of interest here, the configura¬ 
tional entropy Sc[f] is [ 1(1 

Sc[f] = ~ J m\n[f(k)]d d k, ( 10 ) 

where /(k) = /(k)//(k) max and /(k) max is the max- 
imum fraction, in many cases of interest given by the 
zero mode, k = 0, or by the system’s longest physi¬ 
cal mode, | fc m i n | = 7 t/R. This normalization guaran¬ 
tees that /(k) < 1 for all physical values of k. We call 
<r(k) = —/(k) ln[/(k)] the configurational entropy den¬ 
sity. In this paper, we will compute the configuration en¬ 
tropy from the energy density p(r) of the self-gravitating 
object. The choice of the energy density is the most nat¬ 
ural, given that it is a spatially-localized function that 
encapsulates all the relevant physics and boundary con¬ 
ditions describing the stellar configuration. The task at 
hand is thus to solve the relevant Einstein’s equations 
to obtain the equilibrium configurations in terms of p(r) 
and then use p(r) to compute the CE as a function of 
the star’s central density, p(r = 0) = pQ. We start with 
the simplest case, Newtonian polytropes modeling cold 
white dwarfs. 


IV. COLD WHITE DWARFS AND THE 
CHANDRASEKHAR LIMIT 


where the constant K depends on the entropy per nu¬ 
cleon and chemical composition. No heat flow through¬ 
out the object requires 7 to be the adiabatic index, de¬ 
fined as the ratio of the heat capacities of the fluid at 
constant pressure and volume. Small mass, stable non- 
relativistic white dwarfs are well-modeled by 7 = 5/3 

j-2 / o 2 \ 5/3 

and K = 15y ^ — 2 ( J , where m e ^) is the electron 

(nucleon) mass, and p ~ 2 is the number of nucleons per 
electron. The largest mass white dwarfs are modeled by 
/ 2 \ 4 / 3 

7 = 4/3 and K = ( 7777 J ? the well-known Chan¬ 

drasekhar limit nun- The binding energy for polytropes 
with Q nucleons, Eb = M — Qm at, can be written as 
Eb = — 5 w ^ ere M is given by Eq. lil There is 

a clear stability boundary at 7 = 4/3 where Eb changes 
sign. We will show below that the configuration entropy 
captures the same boundary. 


Solutions to Eqs. ITT] and 12 must satisfy p( 0 ) = p 0 


and p'(0) = 0 , and are found introducing new vari¬ 
ables p(r) = po<9(0 V(7-1) an d £ = r / a , with a 2 = 
^GPf- i) Po 7 ~ 2 ^ Equation 11 then becomes the Lane- 
Emden equation with boundary conditions 6(0) = 1 and 
0 '( 0) = 0 , 



2^ + 0l/(7-l) =O . 


(13) 


Solutions were obtained via a 4 t/l -order Runge-Kutta 
method with step size 10 -3 . The CE is computed from 
the energy density using Eq. [lO] Since polytropes have a 
finite radius (where p(R) = 0 or, equivalently, $(£#) = 0 , 
with = ii/a), the /c-integration is in the interval 
k G [fcmin = 7 r/R, 00 ). This ensures that only modes 
with wavelengths smaller than the polytrope contribute 
to the configurational entropy DU- In Fig. [l]we plot the 
normalized modal fraction /(|k|) for sample values of the 
polytropic index 7 . The infrared cutoff is at k m i n = 7 t/R, 
defined by the star’s radius. The configurational entropy 
for the various polytropes is computed using this modal 
fraction to integrate the CE density as described in sec¬ 
tion III. 

We next present a simple scaling argument relating 
the stellar mass to its configurational entropy. Using the 
dimensionless variables 6 (g) and £ in Eq. [T] we obtain, 


Newtonian polytropes are obtained from the hydro¬ 
static equation (setting the general relativistic correc¬ 
tions to zero in Eq. § Eld 


d I" r 2 dp(r) 
dr [p(r) dr 


— 47 t Gr 2 p(r). 


( 11 ) 


Eq. |TT] is supplemented by a general polytropic equation 
of state 


P = Kp\ 


( 12 ) 


r R 

M = 4n p(r)r 2 dr 
Jo 


4irp 0 a 3 f 0 1/(7 
Jo 


3 (3 7 —4)/2 

oc p 0 a oc pQ . 


(14) 


Using the dimensionless variables in the Fourier trans¬ 
form of the energy density, we can express the modal 
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FIG. 1. (Color online.) Normalized modal fraction /(|k|) for 
sample values of the polytropic index 7 . From left to right, 
7 = 1.2,1.4, and 1.7. 


FIG. 2 . (Color online.) Configurational entropy times p^ 1 
(continuous line) and mass (dotted line) versus polytropic in¬ 
dex 7 for po = p c . 


fraction as 


f ( k ) = 


h(ak) h(ak) 

Wf) = H£) 

h(ak) 


(15) 


C{ 7)’ 

where C( 7 ) is independent of po> and h(ak) is 

r^R 

/o 


h(ak) = 


f 0 V (7 !)(£) exp(i<a& • £)£ 2 d£ 

io 


(16) 


The configurational entropy is then, 


5 = —47t 


f 

J kry 


h(ak) 1 ^ h(ak) \ 7 2 


= —47TQf 


-3 


C(7) ‘ 0g 
h{ki) 


C(l)J 


k z dk 


r ««), 1 • 


C( 7 ) 108 (sm) ^ 


(17) 


where k = ak, so that 


7 r/^#. We thus obtain, 


S/00 


oc a 


Vo 1 


ap* 4 " 3 ^ 2 . 


(18) 


Comparing with Eq. 


14 


we see that oc M -1 . Note 
that at 7 = 4/3 the quantity -S '/// 1 is independent of po, 
consistent with a boundary in the star’s stability [ 6 ]: as 
is well-known, stars with 7 < 4/3 are unstable, while 
stars with 7 > 4/3 are stable. 7 = 4/3 defines an in¬ 
stability ridge for the family of stellar configurations, as 
we show next by exploring how both the mass and the 
configurational entropy vary with central density and 7 . 

The mass and configurational entropy Sp$ 1 are shown 
as a function of polytropic index 7 in Fig. i with po = pc, 
where p c is a fiducial value for the critical central density, 


which can be computed for a few specific values of 7 . For 
example, for 7 = 5/3, p c = 0.97 x 10 6 p gm/cm 3 , where 
p ~ 2 is the number of nucleons per electron [6j. 

In Fig. [3] we show the contour plot of the stellar mass 
as a function of po/p c an d 7 , where the existence of a 
saddle ridge at 7 = 4/3 is clear. In Fig. [I] we show the 
contour plot of the quantity Sp^ 1 as a function of po/pc 
and 7 . The reader can verify that the shadings are ap¬ 
proximately reversed for the two plots (there are small 
deviations due to the gamma-dependence of the relevant 
quantitities), illustrating qualitatively the inverse scal¬ 
ing between mass and configurational entropy discussed 
above. 

In Fig. [5] we plot the configurational entropy versus 
polytropic index 7 for various choices of cutoff for fc m i n . 
We do this to illustrate the sensitivity of the results to the 
choice of cutoff and to establish that it is best to choose 
what is physically more natural, that is, no arbitrary fine- 
tuning of cutoff, and thus k m [ n = tt/R. Even if the results 
are not perfectly accurate with this choice, they lie within 
a few percent from the critical points for CE: the maxi¬ 
mum of CE lies at 1.3% from 7 = 4/3- the polytropic 
index for an ultra-relativistic white dwarf - and the min¬ 
imum of CE at 2.1% from 7 = 5 / 3 , the polytropic index 
for the most stable nonrelativistic white dwarf. In par¬ 
ticular, the CE has a maximum for 7 = 1.316. Since the 
mass decreases with 7 for fixed po / p c (cf. Fig. [2|, the CE 
gives an upper bound on the maximum mass for stabil¬ 
ity, M max . This mass approximates the Chandrasekhar 
mass (Mch an d ra ) a t 7 — 4/3 as Tf max — E0373A4ch anc j ra , 
so within 3.73% of the correct value. We also correct a 
typo in Ref. HQ in the vertical axis, which should read 
Sa 3 and not SR 3 . 

We see that the CE accurately captures the stability 
properties of Newtonian polytropes. This gives us con- 
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FIG. 3. Contour plot for the mass of polytropes as a function 
of the central density po/p c and polytropic index 7 . There is 
an instability ridge-a saddle line-for 7 = 4/3. The vertical 
bar specifies the values in units given at its bottom. 


FIG. 5. (Color online.) Configurational entropy versus poly¬ 
tropic index 7 for polytropes. We display results for several 
choices of cutoff for k m in- 



FIG. 4. Contour plot for the configurational entropy of poly¬ 
tropes as a function of the central density po/p c and poly¬ 
tropic index 7 . There is an instability ridge-a saddle line-for 
7 = 4/3. The vertical bar specifies the values in units given 
at its bottom. 


fidence that we can use similar methods to examine the 
stability of general relativistic compact objects, which we 
do next. 


with an Oppenheimer-Volkoff (OV) equation of state, 
where neutrons are treated as a pure ideal gas [T2] . Much 
work has been done in the past decades extending the 
results of OV to more realistic situations, where calcu¬ 
lations try to incorporate a variety of effects taking into 
account the role of the strong nuclear force at the star’s 
core. These approaches are treated in many books and 
reviews, such as those listed in references mmi- Our in¬ 
terest at this point is not to explore different equations of 
state modeling a neutron star interior, but how the effects 
of general relativity, in particular their impact on a star’s 
stability, are reflected in its equivalent configurational en¬ 
tropy. Can we obtain information about a compact rel¬ 
ativistic object’s stability from its information-entropic 
complexity? 

Considering a gas of particles with rest mass /iq obey¬ 
ing Fermi-Dirac statistics, the related equation of state 
may be written in parametric form as fl2] 


p = iC(sinh t — t) 

p = ^iF(sinht — 8sinh ^-t + 3£), 
o Z 

where K = /4h 2 and 


t = 4 log 




(19) 


( 20 ) 


V. NEUTRON STARS WITH 
OPPENHEIMER-VOLKOFF EQUATION OF 
STATE 


As we move into general-relativistic objects, we start 
with a simple but representative model, neutron stars 


with kp being the maximum momentum in the Fermi 
distribution, related to the particle number density n as 
n = kp/37r 2 h 3 . We follow OV and introduce a new mass 
function variable u(r) (equivalent to the mass density 
function GM(r ) defined in Section II), 

1 = 1- 2 u/r, (21) 
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so that u obeys 

du/dr = Airpr 2 . 

Einstein’s equations (Eqs. [5]) become: 

du 9/ . , 

— = r (sinh t — t) 
dr 


dt 

dr 


sinh t — 2 sinh \t 


r (r — 2?x) cosh t — 4 cosh + 3 
^r 3 (sinht — 8 sinh ^-t + 3 1) + u 


( 22 ) 


(23) 


where K — 1 / 47 r, and the units of length a and mass 6 
have been fixed as 


again k m j n = tt/R < k < oc, reflecting the fact that 
neutron stars have well-defined radii where p(R) = 0. 
The results are shown in Fig. [7| where we can see that 
the quantity Sp$ 1 has a minimum (p m in) near the crit¬ 
ical equilibrium value of the central density (p c ) where 
the stellar mass is a maximum. The configurational en¬ 
tropy is multiplied by the inverse central density to have a 
quantity that scales with dimensions of inverse mass.The 
inset shows the results in more detail near the CE min¬ 
imum at p m i n = 0.6196a -3 and thus within 5.3% from 
p c . We can translate the value of p m i n to an estimate of 
the critical mass based on CE, thus establishing a bound 
in the critical OV neutron star mass with accuracy of 
0.58%. (The flatness near the mass function maximum 
helps.) 


1 ( h \ 3/2 c 
^ VMo c) ( Mo G) 1/2 ’ 



(24) 


We can now solve numerically Einstein’s equations 
with boundary conditions 


u(r = 0 ) = 0 ; t(r = 0 ) = to 
u(r = r b ) = u b \ t(r = r b ) = 0 , 


(25) 


so that p(r = r b ) = p(r = r b ) = 0 and r b is the radius of 
the star and u b is its mass. Results are thus parameter¬ 
ized in terms of to, related to the star’s central density po 
by Eq. [l9j In Fig. [ 6 ] we plot the mass of the OV neutron 
star for values of the central density parameter po- Stars 
with p 0 > p c are perturbatively unstable to gravitational 
collapse, as is well-known. 




FIG. 6. OV neutron star mass vs. central density po. Stars 
with po > p c — 0.5886a -3 are known to be unstable to gravi¬ 
tational collapse. 

As with the Newtonian polytropes, we compute the 
configurational entropy using the energy density of the 
equilibrium configurations. The range of integration is 


FIG. 7. Configurational entropy times p^ 1 for the OV neutron 
star vs. central density po- Stars with po > p c = 0.5886a -3 
are perturbatively unstable to gravitational collapse. The in¬ 
set shows the result near the CE minimum. 

The reader should not confuse the results of Section IV 
for Newtonian polytropes, where the prediction for the 
Chandrasekhar mass was given at the maximum of the 
CE with respect to the polytropic index 7 , with the re¬ 
sults here, where the estimate for the critical mass comes 
at the minimum of the CE with respect to central den¬ 
sity po- The OV equation of state is only well-modeled 
by a polytrope in the non-relativistic limit for neutrons, 
with 7 = 5/3. For this value of 7 , the star’s mass is 
a decreasing monotonic function of po 0. The maxi¬ 
mum mass in the mass vs. central density plot appears 
only when general relativistic effects are included. From 
Fig. 0 that the configurational entropy offers a reliable 
measure for the stability of OV neutron stars, provid¬ 
ing an accurate estimate for the critical mass. We also 
verified numerically that neutron stars obey the same ap¬ 
proximate scaling as Newtonian polytropes, in that the 
product Spq 1 M is nearly constant with po- The result 
is shown in Fig. [ 8 j This inverse scaling justifies why the 
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critical mass, being a maximum, correlates well with the investigate whether the configurational entropy of boson 
minimum of the configurational entropy. stars can furnish information about their stability. 



FIG. 8. Behavior of the quantity Sp^M as a function of 
central density po for OV neutron stars. The approximate 
linear scaling persists for a wide range of central densities 
The inset shows the result near p c = 0.5886a -3 , the critical 
value for stability. 

We now extend our approach to another class of gen¬ 
eral relativistic bound objects, boson stars. 


VI. BOSON STARS 

Boson stars are self-gravitating spheres of scalar fields 
Eununi. These hypothetical objects are possible both 
in the Newtonian and general-relativistic limits, and for 
free and self-interacting fields. Due to their remarkable 
properties, boson stars have attracted much interest over 
the past decades [ 8 |. Being made from self-gravitating 
spin -0 bosons, these objects are supported against grav¬ 
itational collapse from Heisenberg’s uncertainty pressure 
and, when applicable, from repulsive interactions among 
the particles m- Indeed, for a free complex scalar 
field with [/(l)-conserving charge Q, these objects have 
masses M (M?i/ m) and radii R ~ 1/m, where m is 
the mass of the particle excitation of the field, and Mp\ is 
the Planck mass mmm- If a repulsive self-interaction 
is added, the mass scales as M ~ A 1 / 2 (Mp 1 /m 2 ) m- 
Furthermore, boson stars have many qualitative similar¬ 
ities with neutron stars, with a maximum mass M mark¬ 
ing the stability boundary against radial perturbations 
na cm. Stars made of real scalar fields are also pos¬ 
sible, but the configurations are time-dependent, known 
as oscillatons [ 22 ] . Given the many similarities, and the 
potential applicability of boson stars to many questions 
of current interest, from being dark matter candidates 
to serving as exploratory tools probing the boundary be¬ 
tween classical and quantum field theory, we will now 


A. Formalism 


For completeness, we briefly review the essential for¬ 
malism to find boson stars. Consider the action 

s = j d 4 xy/=g {—( 26 ) 

where C is the Lagrangian density, 

£ = g^d^dvcj)* - m 2 \<t>\ 2 - ^|</>| 4 . (27) 

We write the spherically-symmetric complex scalar field 
as 0 (r, t) = ^{r)e~ lujt , where <F(r) is real and has no 
nodes. This means that we will only be investigating 
here the stability properties of boson stars in their ground 
state. Stars can be found in excited states and their de¬ 
cay properties have interesting consequences, including 
the generation of gravitational wave bursts [23]. We de¬ 
fine the dimensionless variables x = mr and t = mt. 
Primes are derivatives with respect to x. We also absorb 
the dimensionless frequency Cj = uj/m into the metric co¬ 
efficient B, B = B/uj 2 and define the dimensionless field 
cr(x) = ^(x)/VSttG. (Henceforth we suppress tildes.) It 
proves convenient to rewrite the dimensionless coupling 
constant A as 15] 


A = A 


Mpi 

87r m 2 


(28) 


With these definitions, variation with respect to the met¬ 
ric of Eq. [I] and with respect to the scalar field gives Ein¬ 
stein’s equations (Eqs. [ 2 | and the Klein-Gordon equation 
as 


A' = xA 2 


B' = xAB 


cr' 2 1 2 

-A +{ B +1) ° 


+ 


a' 2 .1 

A + ^B 


l )^ 2 - ^ 


4~ 
Act 4 
4 


- -(A- 

X 

+ -(A- 
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1) 


1) 


"2 

1 , 

( B' 

A'V 


iy 1 ^ 

1 A , 

_ 

T — 



a' — A 


CT-CT 3 

X 

2 ' 


~AJ. 


l\B J 

1 2 


( 29 ) 

These equations are solved for the boundary conditions 
A(0) = 1; B( oo) = 1; cr(0) = <to; cr(oc) = 0; cr'(0) = 0. 
Note that solutions are parametrized by the central value 
of the scalar field cro, which determines the star’s central 
density. 

In Figure [9] we plot the boson star mass and conserved 
charge for the free field case (A = 0). Note that the 
maximum mass is also where the binding energy E^ is 
maximal, where — M — Qm (lower line). As shown 
in Refs. mm the maximum mass is also the stability 
boundary for the boson star. This is also the case for the 
interacting case, A 7 ^ 0 mm- Note also that stars with 
00 > 0.540 have E^ > 0, and are thus unstable to fission. 





















FIG. 9. Boson star mass (continuous line) and conserved 
charge Q (dot-dashed line) vs. central value of the scalar 
field (Jo. We also show the binding energy Eb (dotted line). 
Stars with cro > cr c = 0.271 are known to be unstable to 
gravitational collapse. 


B. Configurational Entropy for Boson Stars 

We now compute the configurational entropy for bo¬ 
son stars from Eq. [lO] using the Fourier transform of the 
energy density as we did with Newtonian poly tropes and 
neutron stars. This means that for each value of cro we 
find the solution of the coupled Einstein-Klein-Gordon 
system of equations and use it to compute the star’s en¬ 
ergy density p(r). We do this for several values of the 
scalar self-coupling. Note that since the scalar field only 
vanishes at spatial infinity, boson stars don’t have a spe¬ 
cific radius where the energy density and pressure vanish. 
We thus don’t use a momentum cutoff, computing the CE 
for all momenta 0 < |k| < oo. The results are shown in 
Figure [lO] as a function of the field’s central value cro for 
different values of the coupling A. The vertical lines de¬ 
note the critical value of the field (<r c ) beyond which the 
star is unstable under radial, charge-conserving pertur¬ 
bations. It is apparent that these lines are very near the 
minima of CE for all values of A. Just as for OV neu¬ 
tron stars, the configurational entropy provides a reliable 
bound on the star’s stability. 

Specifically, we can use the minimum of the configura¬ 
tional entropy to obtain a bound on the star’s maximum 
mass, as was done previously for neutron stars. The re¬ 
sults are summarized in the Table below for different val¬ 
ues of A. The third column is the value of the critical 
mass obtained from the CE, while the last column gives 
the percentual error of the estimate. 


Given the qualitatively similar results between neutron 


A 

M cr it 

Merit 

A(%) 

0 

0.6330 

0.6324 

0.10 

10 

0.7863 

0.7845 

0.25 

50 

1.2450 

1.2367 

0.66 

100 

1.6522 

1.6351 

1.04 


TABLE I. Comparison of maximum stable mass for boson 
stars obtained from traditional perturbation methods (second 
column, Merit) and from the configurational entropy (third 
column, M^f) for different values of the scalar field self¬ 
coupling A. The fourth column list the percentual error of 
the estimate using the CE. 





FIG. 10. The configurational entropy for boson stars mul¬ 
tiplied by inverse central density Sp J 1 as a function of the 
scalar field’s central value cro for different values of the scalar 
field coupling A. The dashed vertical line denotes cr c , the in¬ 
stability boundary for the star under radial perturbations. As 
in the case with neutron stars, the CE provides a reliable es¬ 
timate for the critical mass, with precision better than ~ 1%. 
Values of maximum stellar masses for different values of A are 
listed in Table [I] 


stars and boson stars, we should be able to show that the 
quantity Scr^^M is approximately constant with respect 
to the central field value cro, leading to the inverse scaling 
S<Jq A ~ M -1 . The fourth power of the field ensures the 
same dimensionality as when the energy density is used, 
as was the case with polytropes and OV neutron stars. 
The result is shown in the inset of Fig. El where it can 
be seen that an approximate linear scaling holds in the 
region near and above the critical energy density. The 
main plot shows the second derivative of the quantity 
approaching zero near the critical value cr c = 0.271. In 
Fig. [T2] we present the equivalent results for stars with 
A = 100, where a similar scaling holds. 
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FIG. 11. The inset shows the quantity Sc/)q 4 M as a function 
of the scalar field’s central value cro for A = 0. The main 
plot shows the second derivative of Sc/)q 4 M with respect to 
cro, displaying the same saddle ridge behavior found for poly¬ 
tropes. 



FIG. 12. The inset shows the quantity S(/)q‘ 4 M as a function 
of the scalar field’s central value cro for A = 100. The main 
plot shows the second derivative of S(/)q 4 M with respect to 
cro, displaying the same saddle ridge behavior found for poly¬ 
tropes. 

VII. CONCLUDING REMARKS 

We investigated the stability properties of a variety 
of self-gravitating compact objects using the recently- 


proposed configurational entropy m , a quantity that 
computes the relative weights of different Fourier modes 
making up a given configuration inspired by Shannon’s 
information entropy m- We extended previous results 
for Newtonian polytropes, where the Chandrasekhar 
mass for ultra-relativistic white dwarfs is estimated to 
within 3.73% of the correct value, to fully general- 
relativistic neutron stars modeled with an Oppenheimer- 
Volkoff equation of state and to boson stars made of 
self-interacting complex scalar fields. Using the energy 
density of the configurations to compute their respective 
configurational entropy, we were able to obtain predic¬ 
tions to the critical stable mass with precision better than 
one percent for all these objects. We have further shown 
that an inverse scaling relation holds between the star’s 
configurational entropy and its mass near the critical re¬ 
gion and beyond. This scaling helps clarify why a critical 
value for the mass is reflected in a critical value for the 
configurational entropy, although we are still pursuing a 
first-principles derivation relating the two quantities. 

We are current investigating two related questions. 
First, we are computing the configurational entropy of 
excited states of boson stars [23] in order to relate 
their decay and gravitational radiation emission to their 
configurational-entropic properties. Preliminary results 
indicate that the configurational entropy grows with the 
quantum numbers labeling excited states (n,£, m), as is 
the case with simple quantum mechanical systems. It will 
be interesting to see whether the configurational entropy 
will provide information to resolve energy-degenerate 
states, i.e., states with the same binding energy and dif¬ 
ferent quantum numbers. We are also investigating the 
evolution of the configurational entropy during gravita¬ 
tional collapse. We expect that as the star becomes more 
localized its configurational entropy will increase. An im¬ 
portant question is to determine whether the configura¬ 
tional entropy reaches a maximum when the event hori¬ 
zon forms and whether there is a relation between this 
hypothetical maximum value and Bekenstein’s entropy 
based on the surface area of the black hole [24]. 
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